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ABSTRACT 

The OGLE project led to discovery of earlier unknown forms of multiperiodic pulsation in 
Cepheids. Often, the observed periods may be explained in terms of simultaneous excitation of two 
or rarely three radial modes. However, a secondary variability at about 0.6 of the dominant period, 
detected in a number of the first overtone (10) pulsators inhabiting the Magellanic Clouds, seems 
O j to require a different explanation. After reviewing a possibility of explaining this signal in terms of 

radial and nonradial modes, I find that only unstable modes that may reproduce the observed period 
ratio are f-modes of high angular degrees (£ = 42—50). I discuss in detail the driving effect behind 
the instability and show that it is not the familiar opacity mechanism. Finally, I emphasize the main 
difficulty of this explanation, which requires high intrinsic amplitudes implying large broadening of 
spectral line. 
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1. Introduction 

When a secondary periodicity is detected in a Cepheid, the initial hypothesis 
is always that this is due to excitation of another radial mode. There are some 
theoretical and observational arguments against interpretation in terms of nonradial 
modes. None of the arguments is compelling, especially if the amplitude of the 
secondary variability is low. However, the interpretation in terms of radial mode is 
the most attractive because only then the secondary periodicity yields immediately 
a precise constraint on stellar parameters. The strength of information contained in 
data on two radial mode periods is best seen in the impact of the seminal Petersen's 
(1973) paper on stellar physics. The huge mass discrepancy revealed in the diagram 
Ps/Pl vs. logPi known as the Petersen diagram prompted the major stellar opacity 
revision. Such diagrams are still commonly used to confront data with models. 

In 1973 the list of double-mode Cepheid contained 8 objects. The number of 
such objects has been considerably increased thanks to massive photometry of stars 
in the Magellanic Clouds within the framework of the microlensing projects: MA- 
CHO, OGLE, and EROS. In particular, the OGLE-III Catalog of Cepheids in the 
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LMC lists 61 objects pulsating simultaneously in fundamental and first overtone 
modes (F/10), and 203 in two consecutive overtones (10/20). Similar numbers 
(F/10 - 59, 10/20 - 215) of double mode Cepheids are listed in the catalog for 
the SMC (Soszyriski et al. 2008). The catalog contains also some rare objects with 
frequencies which may be interpreted in terms of different combinations of radial 
modes. Exceptionally valuable are objects where three radial modes are identi- 
fied. Data on the three periods allow for an independent determination of the dis- 
tance to Magellanic Clouds (Moskalik and Dziembowski 2005). With a precisely 
determined distance, massive photometry yields strong constraints on models of 
individual double-mode Cepheids. 

In a number of Cepheids, the period of the secondary variability is incompatible 
with the radial mode interpretation. First such cases were discovered by Moskalik 
and Kolaczkowski (2008, 2009) in their analysis of OGLE-II data for the LMC. 
They found such periods in 42 objects with dominant mode being first overtone 
and named them FO-V Cepheids. In most of the cases the secondary periods were 
close to the dominant modes and were interpreted as an analogue of the Blazkho 
effect in RR Lyr stars. However, in 7 objects, the secondary periods were between 
0.6 and 0.64 of the dominant (first overtone) period and required a different inter- 
pretation. After rejecting the possibility that the secondary periodicity is associated 
with another radial mode, the authors concluded that the high-frequency signal is 
due to a nonradial modes but noted the problem with driving such modes. The 
difficulties of this interpretation were also briefly discussed by Dziembowski and 
Smolec (2009). It is interesting that though the Cepheid masses (1-2 M ) deter- 
mined by Petersen (1973) were by at least by a factor of 3 lower than accepted at 
his time, still interpretations other than simultaneous excitation of F and 10 modes 
were never considered in published papers. 

The subsequent phase of the OGLE project led to the discovery of many more 
Cepheids with the dominant first overtone and the period ratio in the 0.6-0.64 
range. The OGLE-III Catalog provides data on 29 objects of this type in the LMC 
and 139 in the SMC. They seem most puzzling Cepheids posing interesting ques- 
tions. First of all we would like to know what is the nature of the secondary pe- 
riodicity. We also would like to understand why they are so much more frequent 
in the SMC than in the LMC and why there is apparently no counterpart of this 
phenomenon among Cepheids pulsating in the fundamental mode. 

In Section 2 1 review the data and reconsider the two-radial mode interpretation, 
taking into account data on periods and reddening-free magnitudes. On the model 
side, linear instability of considered modes is set as the requirement. Interpretation 
of the secondary periodicity in terms of nonradial mode excitation is studied in 
Section 3. Section 4 is devoted to explanation of driving mechanism acting in 
various modes. Visibility of high-degree modes which reproduce the whole range 
of puzzling frequencies is discussed in Section 5. 
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2. The lO/X Cepheids in the Magellanic Clouds 

2.1. TheLMC 

The difficulty of the two-radial mode interpretation for the 10-X Cepheids in 
the LMC is demonstrated in Fig. 1. The data shown there are from Soszyhski et 
al. (2008) and they include in addition to the two periods, the Wesenheit index, 
which is a reddening-free stellar luminosity measure evaluated from mean / and V 
magnitudes with the standard expression 

W/ = /-1.55(V-7). (1) 

The envelope models and radial pulsation periods were calculated with the 
standard Warsaw codes (see, e.g., Moskalik and Dziembowski 2005) for stellar 
parameters leading to unstable 10 modes in the period range of the 10/X pul- 
sators. Two sets of abundance parameters were adopted (Z,Y) = (0.008,0.256) 
and (0.004,0.24) , which are most often used for modeling young stars in the LMC 
and SMC, respectively. For each selected value of Z, M (mass in solar units), and 
L (luminosity in solar units) two values of the effective temperatures were selected. 
The higher one corresponding to the onset of the 10 instability defines blue edge 
(BE). The one by 0.04 dex lower was adopted as the red edge (RE) of the consid- 
ered ranges of T e ff . With our code we cannot determine the red edge. The choice of 
0.04 is in a crude agreement with results of nonlinear modeling of Cepheid pulsa- 
tion with time-dependent versions of the Mixing Length Theory (MLT, Feuchtinger 
et al. 2000) and the observed range of the V — I color (see Fig. 3). 

The MLT parameter a = 1.5 was adopted for this plot but the value of a has 
only marginal influence on calculated periods. The growth rates are more affected. 
With a = 0.5, BE and RE are shifted by about —0,015 dex. An increase of a 
above 1.5 does extend the instability range of the third overtone but even at a = 2.5 
this modes remains stable in the relevant range of stellar parameters. The ranges of 
luminosity were chosen to cover the whole range of logPio- Models used in this 
plots are realistic in the sense of physics but do not follow from stellar evolution 
calculations. In fact, at the considered Z values no such calculation predicts enter- 
ing the instability strip at M = 3 in the He burning phase. At M = 4 some do but 
only models that include overshooting reach logL > 3 in this phase. 

In the upper panel of Fig. 1, we may see that the objects split into two nearly 
equinumerous groups with period ratios in the [0.599,0.607] and [0.619,0.627] 
ranges. The objects occur in a rather narrow period range, extending from 0.24 
to 0.41 days. One of these objects (shown with the encircled symbol) is an F/IO 
pulsator. The P\o/Pfo = 0.714 ratio at log/V = 0.506 is well within the range 
of F/IO Cepheids in LMC but the A F /Aio = 0.04 is exceptionally low. There is 
nothing that significantly distinguishes this object from the rest in Fig. 1. 
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Fig. 1. Comparison of the periods and the Wj luminosity indices for the 10/X Cepheids in the LMC 
from the OGLE-III catalog with model values. Observational data are shown with symbols. The 
encircled symbol in both panels refers to the object where the fundamental radial mode is excited 
along with the first overtone. The upper panel shows Petersen diagram. The symbols are based on 
the measured P iq and Px periods. The calculated periods refer to the first and third overtones in the 
envelope models with two indicated metal abundance parameters, Z, shown with lines. Individual 
lines cover the range of T e ff of the first overtone instability strip (see text for explanation), at specified 
mass, M, and luminosity, L. The values of M and logL, expressed in solar units, are given above 
each pair of the lines. The lower panel shows the period-luminosity relation employing the Wesenheit 
luminosity index, W/,(seeEq. 1). 



There is no doubt regarding 10 identification of the dominant mode. The only 
radial mode that may be contemplated for association with Px is 30 because for 
realistic stellar models Pjo/Pyo ^ 0.77 and P40/P10 ^ 0-55 in the range of Pio 
shown in Fig 1. At specified P\o the P30/A0 ratio depends both on Z and M . To 
disentangle two effect we need data on luminosity. 

The model values of Wj shown in the bottom panel were calculated with the 
bolometric corrections and color indices are taken from Kurucz (2004) assuming 
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the 18.5 mag distance modulus to LMC. We see in the lower panel that the depen- 
dence on Z is very weak and that models of the same mass lie nearly on a single 
line in the logPio— Wj diagram. Thus, the position in this diagram is a rather 
clean probe of stellar mass. We see that the majority of the 10/X Cepheids in the 
LMC have masses close to 4. Only at short periods they are significantly lower 
but higher than 3. At M ps 4 the plots in the upper panel leave the possibility of 
the 03 interpretation for the secondary mode only for one object, which is near 
logPio = 0.4. Some of the data are reproduced by models with M m 3 but these 
are inconsistent with positions in the lower panel. In many cases the fit requires M 
significantly lower than 3. An independent argument against the 30 identification 
for the secondary mode comes from stability calculations which show that in all 
considered models this mode is stable. Thus, even for the object at log Pi o ~ 0.4 
such identification is excluded. In fact, it is the most distant from the instability 
range. 

2.2. The SMC 

There are nearly five times as many 10/X Cepheids in the SMC than in the 
LMC and, as we may see in Fig. 2, the range of log Pi o where the objects occur 
is by a similar factor wider. The lack of objects with log Pi o < 0.2 in the LMC 
may be explained by much smaller total number of the 10 Cepheids in this galaxy 
(see Fig. 9 of Soszyhski et al. 2010). However, even at longer periods, the number 
of the 10/X is by a factor of 2.4 larger though, beginning with log Pi o = 0.3, the 
number of 10 Cepheid in the LMC is greater than in the SMC. The differences in 
the incidence of the 10/X Cepheids between the two clouds are striking. In this 
paper we focus mainly on the SMC data because much larger population reveals 
the pattern which for the less numerous LMC objects may not be visible. 

Table 1 

Three sequences of 10 Cepheids: mean parameters and standard deviations 



sequence 




log Pi o 


Wi 


V-I 


No 


A 


0.610 ±0.003 


0.106 ±0.071 


15.54±0.27 


0.56 ±0.06 


65 


B 


0.623 ±0.002 


0.297 ±0.066 


14.90±0.30 


0.61 ±0.05 


50 


C 


0.638 ±0.003 


0.443 ±0.065 


14.37 ±0.25 


0.65 ±0.04 


24 



There are three well-detached sequences of 10/X Cepheids concentrated around 
Px/Pio = 0.610, 0.623, and 0.638 marked in Fig. 2 as A, B, and C, respectively. 
These ratios are correlated with periods but there is some overlap. Within each of 
the sequence we see a slow decline of the period ratio with Pio . The mean param- 
eters of stars forming the three sequences are listed in Table 1. Only the objects in 
sequence B have their counterparts in the LMC. The objects in sequence A, except 
of one, occur at much shorter periods than the LMC objects with similar period 
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Fig. 2. Comparison of the periods and the Wj luminosity indexes for the OGLE-III data (Soszyhski 
et al. 2010) for the 10/X Cepheids in the SMC with model values. Models used in this plot were 
calculated assuming Z = 0.004, and the M—L relation given in Eq. (2). Each of the six lines extends 
over the 01 instability strip (as defined in the text) at fixed L and M . In the part shown in the solid 
style also 03 modes are unstable. The adopted range of logL, is [2.7,3.45] and the resulting range of 
M is [3.2,5.17]. The rectangle in the upper panel shows the range of the 10/X Cepheids in the LMC. 

ratios. There are no 10/X Cepheids in the OGLE-III data for the LMC with the 
period ratio as high as found in sequence C. It should be stressed, however, that 
Moskalik and Kolaczkowski (2009) in their analysis of the OGLE-II data detected 
three such objects. 

The color-magnitude and period-magnitude diagrams depicted in Fig. 3 show 
that the 10/X pulsators occur over nearly whole range of parameters of the first 
overtone Cepheids in the SMC except that they are absent at shortest periods. There 
are 320 first overtone Cepheids at log Pi o < 0.056 but in none of them the puzzling 
high frequency signal has been detected although the incidence of the 10/X pul- 
sators in the whole population is about 8.5%. This, in part, may be blamed to 
difficulty of detecting this low amplitude signal in fainter objects but is unlikely to 
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account for the whole effect. 
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Fig. 3. The 10 Cepheids in SMC from OGLE-III catalog in the CM and PL diagram. The 10/X 
objects from the three sequences are marked with different symbols. 



Models used for comparison with the data in Fig. 2 satisfy the mass-luminosity 
relation 



logL = 3.05 + 3.6(logM - 0.602), 



(2) 



which is based on Girardi et al. (2000) evolutionary models of stars with masses 
M = 4 and 5 at Z = 0.004. The tracks with M = 3 do not reach the instability 
strip in the helium burning phase. BaSTI (Pietrinferni et al. 2006) models also at 
M = 4 and 5 do not rich the instability strip but the maximum values of L in this 
evolutionary phase is in a good agreement with the relation given by Eq. (2). Sim- 
ilar conclusion was reached by Smolec (private communication), who used public 
domain code MESA (Paxton et al. 2011). Apparently, new evolutionary codes do 
not describe well the Cepheid phase. In the lower panel of Fig. 2, we may see that 
with the adopted mass-luminosity relation the observational data in the PIO— Wj 
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plane are well reproduced. Interestingly, I found it also for the LMC data. Mod- 
els calculated with Z = 0.008, which have lower M at the same L, produce lines 
running well above observational data. This suggests that the occurrence of 10/X 
pulsation requires low metal abundance and perhaps in this way may explain the 
large difference in the incidence of such a pulsation form between the two systems. 

2.3. Problems with Two Radial Modes Hypothesis 

The problem with reconciling the secondary period with P30 for the SMC 10/X 
Cepheid is different at short and long periods. In the first case, it is the large discrep- 
ancy in the period ratios and in the second it is the stability of the third overtone. It 
is not easy to find modification of stellar structure resulting in the required change 
in the period ratios. Any admissible changes in Z cannot do it. A factor of ~ 10 
increase suffices to get agreement for the sequence B but not for sequence A. The 
plots in Fig. 3 do no reveal any anomaly that would suggest nonstandard models for 
the 10/X Cepheids. The absence of these objects in the short period range, which 
presumably is populated by stars crossing the instability strip, may suggest that 
the occurrence of the secondary periodicity is related to some property of the deep 
interior structure. One conceivable feature affecting mainly frequencies of higher 
overtones could be a composition jump at fractional mass above M r /M k, 0.6 but 
it is not expected in the relevant stellar models. 

The stability problem may seem less severe because the determination of the 
instability ranges suffers from the poor understanding of interaction between con- 
vection and pulsation. However, especially for sequence C objects, we are far from 
the period range where the third overtone is unstable. This is seen not only in re- 
sults of calculations but also in the observational data. There is one object in the 
SMC at logPio ~ —0.65, that is, well below Fig. 2 range where almost certainly 
the third overtone is excited together with the first two overtones (see Fig. 4 in 
Soszyhski et al. 2010). In the LMC there are three 01/02/03 Cepheids all located 
near log Pi o = —0.25 (see Fig. 2 in Soszyhski et al. 2008). Again this is far from 
the 10/X range. 

Perhaps the correct interpretation of the 10/X Cepheids requires going beyond 
the linear pulsation theory. However, since the period ratios are significantly differ- 
ent from the ratios of low integers, the only possibility seems to be a modification 
of the mean stellar structure by the 10 presence leading to a modification of 30 
properties. This possibility remains to be studied. 

3. Nonradial Modes and the Puzzling Frequencies 

The problem with nonradial mode excitation in giants is the strong damping in 
the radiative interior which is a consequence of huge values of the Brunt- Vaisala 
frequency, hence very large radial wave numbers. However as Osaki (1977) first 
showed there is a class of nonradial modes that despite strong damping in the core 
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are unstable owing to the same driving effect as radial modes. The instability occurs 
if the modes are sufficiently well trapped in the envelope and this is possible if the 
angular degree is sufficiently high. In the 7 M Cepheid model, the instability 
started at t = 4 and the period close to 10 and at t = 6 and the period close to 
the fundamental radial mode. Frequencies and growth rates were determined for 
envelope models alone upon assuming the running wave inner boundary condition. 
At the same time, I found (Dziembowski 1977) similar unstable modes in RR Lyr 
models. Formal derivation of the inner boundary is given in that paper. There 
is one-to-one correspondence between trapped nonradial and radial modes, thus 
following earlier habits, we will use names f- and p„ -modes. In particular, the 
10 mode corresponds to pi nonradial modes, all having one node in the acoustic 
cavity. 

Thirty years later, Mulet-Marquis et al. (2007) revisited the problem. They con- 
ducted an extensive survey relying on complete stellar models and their advanced 
tools for solving equations for nonradial oscillations. All numerical results pre- 
sented in the present paper were obtained with the same code I developed over 35 
years ago. A comparison with modern results shows that this is adequate. We all 
share the same major uncertainty, that is, the treatment of convection. 

Fig. 5 depicts the relative driving rates 



is the local contribution to the work integral, for selected models of Cepheid en- 
velopes. The value of r\ varies between — 1 (no driving zone) and and 1 (no damp- 
ing zone). This parameter, which was introduced by Stellingwerf (1978), is abetter 
measure of the instability strength than the usual growth rates and also better pre- 
dictors for the nonlinear development. 

Only in the cooler (RE) model at the lowest L do we find unstable low degree 
modes. Unstable are P3 -modes at I = 3, 4, 5. At t = 5 the period is significantly 
shorter than P30 and the period ratio is not far from that found in sequence A ob- 
jects but the instability is marginal. At lower degrees periods are very close to Pio- 
In any of the models considered by us, not only those used in Fig. 4, the p-mode 
instability range does not extend much beyond P30 if this mode is unstable. The 
f-modes are clearly different. They are strongly unstable. The maximum of r\ oc- 
curs in the observed range of Px/P\o but there is a difficulty following from high 
angular degrees, which in this range are between t = 40 and 50. Thus, large pul- 
sation amplitudes are needed to offset effect of cancellation in the mean variability. 
We will return to this difficulty in Section 6. Let us first try to explain why f-mode 
driving is so efficient. 



/ d 3 jw 
f d 3 x\w 



(3) 



where 
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Fig. 4. The relative driving rate for trapped acoustic modes and high-degree f-modes in selected 
models (M = 3.2 at logL = 2.7 and 3.87 at logL = 3.0). In the logL = 2.7 BE (the 10 blue edge) 
model the /-ranges of unstable pi and p2 modes are [6,21] and [4,10], respectively. In the cooler 
model (RE) at the same luminosity the corresponding /-ranges are [5,22], and [3,12]. In this model 
there are unstable p3 -modes at I = 3, 4, and 5. In the logL = 3 BE, model the ranges are [7,15], 
and [5,6] while in the logL = 3 RE model, the ranges are [4,19], and [4,9]. In all models the f-mode 
shown have £ ranging from 16 to 54 but the instability extends well beyond this range. 



4. The Driving Mechanism behind the f-Mode Instability 



That high-degree modes are unstable over wide range of stellar parameters has 
been known since long time ago (Dziembowski 1977, Shibahshi and Osaki 1981). 
The authors of both papers noted the instability extends beyond the blue edge for 
radial pulsation and that the hydrogen ionization is the main site of the driving 
effect but they failed to explain the driving mechanism. My aim in this section is to 
supply the explanation. 

Let us focus on two models with the same parameters at the surface but differing 
in the efficiency of the convective transport. Plots in Fig. 5 and 6 refer to models 
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Fig. 5. Plotted against the local temperature in outer layers are various quantities relevant for under- 
standing driving and damping effects working in the three selected modes. The model is characterized 
by the global parameters M = 4.69, logL = 3.3,and log T e ff = 3.806. It was calculated assuming the 
mixing length parameter a = 1.5. The selected modes are lO and 30 radial modes and the I = 42 
f-mode (f42). Their periods are listed in the left bottom panel. The top left panel shows the opacity 
derivative K7- = (9logK/3log T)p and the thermal time scale in days, % of the envelope above each 
point. The upper right panel shows the actual temperature in the envelope dlnT/ dlnP = V, the 
ratio of actual to radiative gradient, V/V lat j, and the adiabatic gradient V a( j. The middle and pan- 
els on the left-hand side show the absolute values of the radial eigenfunctions for the pressure and 
temperature perturbation, respectively. The bottom right panel shows the phase difference between 
the two eigenfunctions. Plotted in the middle right panel Q (see Eq. 4) is the local contribution to 
driving per unit of the independent variable. 



calculated with a = 1.5 and 0.5, respectively. The choice of a has a significant 
effect on nonadiabatic mode properties. In particular, the relative driving rates for 
the 01, 03, and £ = 42 f-mode (f42) are, respectively 0.29, -0.65, and 0.87 at 
a = 1.5, while at a = 0.5 the corresponding rates are 0.04, —0.99, and 0.93. 
However, the huge difference in r\ , between 03 and f42 is preserved and this is 
true for all models considered in this paper. This may look like a paradox. These 



Vol. 62 



11 



O 





i 1 i i i 1 i i i 1 i i i 1 i 


i i | unxjjt 


10 


r a = 0.5 




5 


r logT 






Km 







~ * ^* ■ 








/ 


-5 


1 i 1 1 1 1 1 1 1 1 1 1 1 


A i 1 i i i 1 i 



1.2 
| 0.8 

CD 

£ 0.6 

CO 

\b 4 
0.2 





1 


II | 1 1 1 | 1 1 1 j 1 II 


1 ■ 1 _ 




±._ 


r~\~ 




1 \ 

- _ac i 
-ac/ra i ; 
ad / ■ 


i i : 

j 1 ; 




i 

i \ 






-~ \_/\ 


■ / \-= 


: | 1 


i i l i i i l i i i l i i i 


i i i i i n 



50 
40 
a 30 
20 
10 


400 



I i i i | i i i | ij_ i | i i i | i i i 1 1_ 



I j I I I | I I I | I I I | I IE I | H 

f42 1.85dL- \ 

300 |- 10 2.90d t Hi - 
30 1.83d ? ) 



200 
100 




I lj.j.Lj.Jju 



If 


i\ 


If 


i\ 


/f 


A 


A 


A 


i 


\ 


t i i 







20 


Of 


10 









-10 




1.5 




1 


a. 




•s- 

1 


0.5 






■s- 







-0.5 



I I I I I I I I I I I I I I I I I I I I I I I 



/ v ' 

i'^'i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




4.8 4.6 4.4 4.2 4 3.8 
logT 



4.8 4.6 4.4 4.2 4 3.8 
logT 



Fig. 6. Same as Fig. 5 but for models calculated with a = 0.5 instead of 1.5. 

two modes have nearly the same periods; there is an additional energy loss in f42 
resulting from the horizontal radiative flux and yet this mode is driven in nearly the 
whole outer part of the star where it is trapped. In Figs. 5 and 6 the driving layers 
may be recognized as places where Q > . This quantity is proportional to w (see 
Eq. 3) and it is given by 



Q = Qn, 



vv V 9iogr 



(4) 



where Qn is the normalization constant, which is the same for all modes in Figs. 5 
and 6, V = — dlogP/ dlogr, while yp{r) and yr{r) are radial eigenfunctions cor- 
responding to hP/P and 8r/r. The angular and time dependence of perturbations 
is assumed in the form F"'exp(— /oof) . The standard normalization y r = 1 at the 
surface was adopted for the radial eigenfunctions. In Figs. 5 and 6 in addition to Q, 
plotted are \yp\, \yr\, and the phase difference ((ly — (})/>, whose sign sets the sign 
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of Q. 

The classical opacity mechanism acts only in the Hell ionization zone and only 
in 10. Note rising Kp between logT = 4.7 and 4.55 and the associated bump of Q 
for this mode. For 03 this is the damping zone because its amplitude varies rapidly 
and the radiative losses dominate. Convection is present in this zone but carries 
a negligible fraction of energy, even at a = 1 .5 . In the layers above the thermal 
time-scale is much shorter than periods of the modes. Thus pulsation becomes very 
nonadiabatic. The common feature of the modes is steep peak of \yp | centered at 
logT ?a 4.05. Its height depends on mode and a but neither on its localization nor 
its shape. In certain range around its maximum, \yp | is roughly determined by the 
equation 

-^-=23(K T -4)y T + C (5) 

diogr 

where C, which includes all remaining terms from the expression for perturbed 
radial component of the radiative flux, is regarded constant. Associated with the 
\yr \ peak is the peak of Q at the same location, except for 30 at a = 0.5 when we 
see there a deep. Always higher a enhances driving effect in radial overtones and it 
is opposite to f-modes, where higher a results in some reduction of T] . Apparently, 
in the latter case the effect of increased \yp | dominates over the effect of increased 
X . Regardless of the value of a the phase difference does not change much in the 
layer around the \yj\ maximum. It growth takes place below, between logT = 4.6 
and 4.2. Thus, the driving mechanism actually works there and not in layers around 
maximum of Q . 

In order to identify the effect causing instability of the high-degree f-modes, we 
transform the expression for Q (Eq. 4) using 



3| '^J.\ =Z (^1) = 1* 



P T ) \ P c P I c P 



5P* 

(OP y ' 



where F denotes the total energy flux per unit area. Then, we separate Q into the 
part arising from the vertical (Qv) and horizontal (Qh) part of 8(VF) . Adopting 
the convective flux freezing approximation in the form 8(r 2 F conv ) = 0, we obtain 

yp is the radial eigenfunction corresponding to 5 (47tr 2 i 7 ra d ) /Z. , and 

*~ ^^(S<— »)• 

The two contributions and the net value of Q in the f42 mode are plotted in Fig. 7, 
where we may see that below logT = 4.2 the horizontal energy losses lead to 
pulsation energy gain {Qh > 0). In order to understand why is it so let us note 
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that {yj + y r W) is the radial eigenfunction for the relative Eulerian perturbation 
of temperature, T'/T . If CO is regarded real, which is not far from true, we have 
Qh x 9^(8/ 5 *r / ). Thus, driving takes place if cooling, which is proportional to 
T , occurs in the low pressure phase and heating occurs in the high pressure phase, 
just like in the well-known case of semi-convective instability in the presence of the 
mean molecular weight gradient. Long-time ago Souffrin and Spiegel (1967) found 
that the same effect may cause instability of the gravity waves also in chemically 
homogeneous layers. In the present situation, Qh > because in the considered 
range of log T , the advection term y r VV dominates in T and the phase difference 
\[f r — \|/ P k, %/2 is largely fixed by the inner boundary condition. 
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Fig. 7. Values of Qy and Qh in models used in Figs. 1 and 2. 



The pulsation energy gain from Qh > is nearly canceled out by Qy < 
because we are in the region where T <C Px and thus there is a tendency toward 
the temporal thermal equilibrium. This changes near logT = 4.2, where the steep 
rise of opacity inhibits the horizontal diffusion and creates the yj bump, which 
enhances the vertical diffusion. Most of the contribution to driving arises there but 
the essential phase lag \\fj — is created below. 



5. Visibility of High Degree Modes 



We have noted in Section 3 that the f-modes that could explain the Petersen 
diagram for the 10/X Cepheid must have rather high degrees. The plots in the 
upper panel of Fig. 8 suggest the I = 42, 46, and 52 are excited in sequence A, B, 
and C, respectively. The fact that only even degree modes are seen is consistent with 
expectations because cancellation causes more rapid decline of the disc-averaged 
amplitudes for odd-degree modes but why multiples of four are preferred is difficult 
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to explain. The p i modes in the same frequency range have £'s by about a factor 
of 2 lower, hence they might appear a better option for the identification. However, 
the ^-values best fitting for the three sequences are odd (19, 21, and 23) hence are 
disfavored by the visibility argument. The arguments discussed already in Section 4 
also speaks against the p i option. In the RE models, the unstable p i at £ = 19 occur 
only for logP 10 < 0.3 , at £ = 21 for logPio < 0. 1 , and at £ = 23 for log^io < 0, 
hence to the left of the great majority of the data point. In the BE models the 
instability ranges are still narrower. They are also narrower in models calculated 
with lower a . 
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Fig. 8. In the upper panel, the Px/P\o period ratio for the 10-X Cepheids are compared with the 
evaluated ratios of P(/P\o at even £'s for the same models as used in Fig. 2. The labels showing 
mode degrees are put between the lines corresponding to the red and blue edges of the instability 
strip. In the lower panel, the observed amplitude ratios in the /-band are compared with calculated 
amplitudes assuming the same intrinsic amplitudes for the f-modes and the lO mode. 



The calculated amplitude ratios shown in the lower panel were obtained assum- 
ing the same rms amplitudes of the photospheric radius variations, 57? rms , for all 
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modes. The numbers were obtained with the use of parameters from Kurucz (2004) 
models and Claret (2000) description of the limb-darkening law. The disc-averaged 
fluxes were evaluated with help of the recurrence relation following Townsend 
(2003). The observable amplitudes of the f-modes decrease slowly with £ because 
in this range the dominant cause of the light changes is the geometrical distortion. 
With the parabolic limb-darkening law, the amplitude is proportional to £~ - 5 if £ 
is even. The decline is also slow with the Claret version. Much larger ratios for the 
RE than BE models follow primarily from much larger surface values of yp for the 
10 mode. 

The comparison with observations suggests that amplitudes of the f- 

modes would have to be of the same order as that of the 10 mode to explain the 
10/X Cepheids. This is a very crude conclusion because it is based on the linear re- 
lation between the flux and radius amplitudes whereas at such intrinsic amplitudes 
we are then well beyond validity of the linear approximation. We should also keep 
in mind that calculated values refer to aspect-averaged amplitudes while the actual 
amplitudes depend on aspect angle, i, and on the mode azimuthal order, m, which 
we do not know. Certainly, the high intrinsic amplitudes of the nonradial modes 
required for this interpretation is a serious difficulty. However, with interpretation 
in terms of the pi modes, we get the amplitude ratios lower only by some 50%. 
Taking into account the other difficulties of this hypothesis discussed earlier, we 
cannot regard it as a viable alternative. 

The difficulty of the interpretation in terms of the high-£ f-modes is best re- 
vealed with the estimate of the radial velocity variations within stellar photosphere, 



with A v = 2ooo7? rms . This expression uses CO 2 = £GMR~ 3 , which is approximately 
valid for these modes. The value of 57? rms may be expressed in terms of the aspect- 
average light amplitude. Keeping only the contribution from geometrical distortion 
of the stellar surface, we get 



where the photospheric radius, R, is expressed in the solar units and Pio in days, 
while h(jj) describes the limb-darkening law. This quantity, evaluated for the same 
models as used in Fig. 2, is shown in the lower panel of Fig. 9. To estimate the 
local radial velocity we may use data on light amplitude, which are plotted in the 
upper panel. The typical velocity amplitudes range from 35 km/s for sequence C to 
70 km/s for sequence A. These are very large numbers. At so high £ 's, the expected 
signature is a nearly constant spectral line broadening. Such large effect could not 
be overlooked. Unusually broad line profiles (AX/X > 10~ 4 ) were observed in 
certain galactic Cepheids. Kovtyukh et al. (2003), report such a phenomenon in 





-i -l 



[km/s/mmag] 



(8) 




four objects, of which three are certainly and one is possibly first overtone pulsator. 
The authors noted significant changes in the line profiles during pulsation cycle 
including moving bumps. Their cautious interpretation is that low degree nonra- 
dial modes are excited in these stars. Indeed, the observed patterns suggest rather 
low £'s. 

The high degree mode instability is found over the whole luminosity range 
of the classical instability strip and its blueward extension (Dziembowski 1977, 
Shibahashi and Osaki 1981). In all the cases the dominant contribution to excita- 
tion arises in the hydrogen ionization zone. The direct evidence that such modes 
are indeed seen in stars comes from observations of the line profile variations in 
8 Set stars. In particular, Kennelly et al. (1998) report a detection of modes with 
i up to 20 in x Peg. Photometry alone may yield only circumstantial evidences. 
For instance peaks at frequencies above the instability range for p-modes could 
be interpreted as being due to high-^ f-modes if they cannot be accounted for in 
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terms of nonlinear effects of lower frequency modes. This was the argument made 
by Daszyriska-Daszkiewicz et al. (2006) regarding the high-frequency peaks in the 
8 Set star FG Vir. However, we can never be sure that the inventory of the lower 
frequency peaks is complete. 

6. Summary and Discussion 

The secondary periodicity in the first overtone Cepheids at about 0.6 of the 
dominant period presents a challenge to stellar pulsation theory. The phenomenon 
is not rare. Soszyhski et al. (2010) found it in 8.5% of the SMC 10 Cepheids. 
The incidence in the LMC is much lower, but still is seen nearly 2.4% of the ob- 
jects, also the range of period is much narrower. Objects exhibiting such secondary 
periodicity are called in this paper the 10/X Cepheids. In the SMC, they form 
three detached sequences centered around period ratios 0.61, 0.62, and 0.64 in the 
Petersen diagram. Except for avoiding the shortest periods, there is nothing that 
distinguishes these objects from the rest of the first overtone Cepheids. 

One F/IO/X Cepheid was found in the LMC but the F-mode amplitude is ex- 
ceptionally low. No object that could be regarded as the F/X Cepheid has yet been 
found. The phenomenon seems restricted to the first overtone pulsators but not to 
Cepheids. Olech and Moskalik (2009) found the secondary periodicity in seven 
RRc stars in go Cen with the period ratios between 0.608 and 0.622. Period ratios 
in the [0.612, 0.632] range were found in three RRc stars observed with Kepler 
space telescope (Moskalik et al. 2012). Earlier, a peak near O.6P10 was found in 
the RRd star AQ Leo by Gruberbauer et al. (2007). 

Parameters of all 10/X Cepheids in the Magellanic Clouds are well constrained 
by Wesenheit index. Results of linear nonadiabatic calculations for envelope mod- 
els consistent with observed parameters seem to exclude interpretation of the sec- 
ondary periodicity in terms of a higher radial overtone. The closest period is that 
of the third overtone but, in most cases, it is far too long. In remaining cases the 
third overtone is stable. It is unstable only at the shortest 10 periods, where the 
period ratio discrepancy is the greatest. It is very difficult to imagine modification 
in stellar models that could save the interpretation of the 10/X Cepheids as double 
radial mode pulsators. 

A possibility that the secondary periodicity is due to a nonradial mode trapped 
in the envelope was considered in Section 3. The conclusion was that only high- 
degree f-modes remain unstable at the puzzling high frequencies. The instability of 
such modes over wide ranges of frequencies and stellar parameters has been known 
for years but has never been satisfactorily explained. In Section 4, we could see 
that the crucial role in the instability is played by the horizontal radiative transport, 
which results in energy gain (loss) of the gas element in the high (low) temperature 
phase. Nothing is known about nonlinear development of the instability 

The three sequences of the SMC 10/X Cepheids corresponding to Px/Pio ~ 
0.64, 0.62, and 0.61 may be associated with single angular degrees I = 42, 46, and 
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50, respectively. It is not surprising that the degrees are even because amplitude 
reduction caused by cancellation effect is stronger at odd £'s. No explanation was 
given why the selected values are multiples of four. In fact, the greatest difficulty 
of this interpretation are large values of £'s. Typical amplitudes of variability in the 
/-band at the puzzling frequencies are in the 4-10 mmag range. To reach such an 
amplitude against the cancellation effect we must postulate the rms- values of 5R/R 
similar to that of the radial first overtone. Excitation of a high-^ mode with so large 
amplitude would manifests itself in a nearly constant, large (35-70 km/s) spectral 
line broadening. No spectroscopic data are available for any Cepheid known as a 
10/X object. However, overlooking the low amplitude signal is more likely than 
the predicted large line broadening. Still, even if the explanation in terms of f- 
modes looks unrealistic, spectroscopic observations of the 10/X Cepheids are most 
desired because data on line profiles may yield a clue to the correct explanation. 
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